High field fractional quantum Hall effect in optical lattices 
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We consider interacting bosonic atoms in an optical lattice subject to a large simulated magnetic 
field. We develop a model similar to a bilayer fractional quantum Hall system valid near simple 
rational numbers of magnetic flux quanta per lattice cell. Then we calculate its ground state, 
magnetic lengths, fractional fillings, and find unexpected sign changes in the Hall current. Finally 
we study methods for detecting these novel features via shot noise and Hall current measurements. 
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The achievement of strongly correlated atomic systems 
[l| which are accurately described by simple condensed 
matter physics (CMP) models |2| has opened many new 
perspectives for research on strongly correlated systems. 
Such atomic ensembles share many properties with con- 
ventional CMP systems. Several of their features, how- 
ever, distinguish them from CMP setups like e.g. their 
extremely long coherence times and the ability to adjust 
Hamiltonian parameters over a wide range and in times 
short compared to the coherence time. These properties 
enable experiments of controlled coherent quantum dy- 
namics |l|, yj , and the exploration of otherwise unacces- 
sible parameter regimes. Due to their highly non-linear 
behavior these studies may lead to striking new features. 

In this paper we consider such a novel situation aris- 
ing if interacting ultracold bosonic atoms in a two di- 
mensional optical lattice are subjected to a simulated 
huge magnetic field. A variety of proposals for realiz- 
ing magnetic fields in optical lattices exist. One can 
e.g. rotate the lattice as suggested in J4J, yj, use laser- 
induced hopping as proposed in 6], or employ an oscil- 
lating quadrupole potential as shown in |jj . All of these 
methods allow the creation of arbitrarily large fields with 
flux quanta per lattice cell a between < a < 1 and 
can be realized with current or near future experimental 
technology. Optical lattice setups are virtually free of 
imperfections and can be cooled to almost zero tempera- 
ture [l| . Furthermore, their large degree of flexibility Q 
allows for the introduction of well controlled amounts of 
disorder [9|. This is in contrast to most CMP systems 
where only values a < 1 are believed to be accessible 
and disorder is difficult to determine and control. Also, 
conventional CMP systems have fermions with a screened 
1/r Coulomb interaction while we consider bosonic atoms 
interacting via a contact potential. Finally, dynamically 
adjustable smooth trapping potentials are used to confine 
atoms in an optical lattice to a certain region in space ja| 
while CMP systems are often bounded by sharp edges. 

For a single particle in a lattice placed into a huge mag- 
netic field with a ~ 1 a fractal energy band structure 
was predicted by Hofstadter |Hj. However, this regime 
was not thought to be experimentally accessible until re- 



cently and is thus not explored in detail. Many studies 
on interaction effects in magnetic fields centered around 
the fractional quantum Hall (FQH) effect for a <C 1 [llj 
where it is essentially unaffected by the presence of a 
lattice. As a ~ 1 the discrete periodic structure of the 
lattice becomes important and we will find novel effects 
arising from the presence of the lattice. We will show 
that the system is described by a bilayer FQH like model 
near simple rational values a c = l/n where n is a small 
integer. The ground state of this system is an FQH state 
with a mean filling factor v = n/(n + 1) and a modified 
magnetic length. We will find striking features of the 
associated Hall current which changes its sign at a = a c . 
We consider ultracold bosonic atoms in the lowest mo- 
tional band of a 3D optical lattice with lattice period 
d. The motion of the atoms along the z-direction is sup- 
pressed by high optical potential barriers [a] . The system 
is thus decoupled into 2D planes and we assume that a 
large effective magnetic field along the z-axis is simu- 
lated by one of the methods introduced in (J, |y, L3| . The 
Hamiltonian in one of the planes is given by 

H = -J^2(e 2maq al q a p - ltq + al q a Ptq - l +h.c.) 

+V(p, q)al q a p . q + -^al q a^ q a p , q a p , q . (1) 

Here a P q (a Piq ) are bosonic creation (annihilation) oper- 
ators, fulfilling standard bosonic commutation relations, 
for particles at a lattice site labeled by its position in 
the xy-plane (p, q) . The onsite interaction strength is U 
and J is the hopping amplitude. The parameter a is a 
measure for the strength of the simulated magnetic field 
[gj . The phase shift from hopping around one lattice cell 
is 27ra, and since a is only defined mod 1 we restrict 
— 1/2 < a < 1/2. Finally, V(p, q) is a trapping poten- 
tial varying slowly on the length scale d. We will mostly 
deal with linear geometries V(p, q) = V{q), where the x- 
dependence of the atomic wavefunctions are plane waves, 
and thus work in Landau gauge. All parameters in the 
Hamiltonian can dynamically be varied via lasers [8j. 

We first discuss the limit \a\ <C 1, in which the length- 
scale of the wavefunctions is much larger than d and we 



can hence use a continuum approximation. We approxi- 
mate a one particle state \ip) = ^ tp(p, q)at q \v&c) with 
|vac) the vacuum state by a continuous wave function 
4>{pd, qd) = ip(p, q)/d 2 . The Hamiltonian acting on <j> is 
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(2) 
where we have defined the effective mass m = l/2Jd 2 
and the cyclotron frequency fl = ira/md 2 . The last term 
has no dynamical effect and is left out in the following. 
For interacting particles at density p, and interparticle 
spacing p^ 1 ' 2 3> d, the continuum approximation is 



ffwE^'i 



J2 s ( x i-x j )S(y l -y j ), (3) 
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where (xi,yi) are the coordinates of particle i and u — 
Ud 2 . By analogy with the solid state FQH we define the 
filling factor v = pw/mil. The 2D Hamiltonian Eq. 
has been studied e.g. in the context of rotating ultracold 
atomic gases for the case of a weak trap 0, Q, u2, QjJ and 
it is well established that for an isotropic trap V(x, y) — 
muj 2 (x 2 + y 2 )/2 the v = 1/2 Laughlin state with a mag- 
netic length t given by £ c = (m 2 fl 2 + m 2 ^ 2 ) -1 / 4 is an 
exact eigenstate. Its energy is N(n 2 + uj 2 ) 1 ^ 2 /2 + N(N ~ 
1)AE, where AE = (ft 2 + lj 2 ) 1 / 2 - fi sw u 2 /2fi and N is 
the number of particles. It is the ground state of the sys- 
tem for u <C ft, umft. Since this state already has exactly 
zero interaction energy, the rest of the bosonic Laughlin 
sequence (y = 1/4, 1/6, . . .) does not occur. It is also 
possible to construct non-circular Laughlin-like states for 
other trap geometries. In general these are only exact 
eigenstates for zero trap strength, but as they are always 
lowest Landau level (LLL) and zero interaction energy, 
they are reasonable trial states for weak traps. For in- 
stance, in linear geometry with a trap V(y) = muj 2 y 2 /2 
the variational energy is minimised for a Laughlin like 
state with magnetic length £i = (m 2 ft 2 + m 2 w 2 /4) -1 / 4 . 

These Laughlin states cannot remain the ground states 
if N gets too large as for N > u/4tt£ 2 AE the trap po- 
tential energy 2NAE to add a particle to the edge of the 
Laughlin state exceeds the interaction energy ~ u/2ir£ 2 
required to add it in the centre of the trap. For larger N, 
denser states are more favourable, and it is believed |lj] 
that a series of incompressible states with half-integer fill- 
ing v = 1, 3/2, 2, . . ., well approximated by symmetrized 
products of Laughlin states called Read-Rezayi states 
[15|. forms. For v > 6 a transition to a compressible 
vortex lattice is predicted [14j . 

We now turn to the case of larger values of a which 
can be reached in optical lattice setups |6J. We re- 
strict our considerations to values of a close to sim- 
ple rationals a c = l/n, where I, n are small integers. 
For a single particle we find from numerical calcula- 
tion shown in Fig. ^ that when a ~ a c the sin- 




FIG. 1: Wavefunction for (a) a = 1/2 and (b) a = 1/3, 
in linear geometry with trap V(y) — muj 2 y 2 /2 and uj — 
10~ 3 7r/md 2 . The dots show the numerical wavefunction and 
the solid curves the analytic approximations for y/d = mod 
n (upper curve) and y/d ^ mod n (lower curve). 



gle particle ground state wavefunction has an approx- 
imate n-site periodicity, suggesting the representation 
ip(np + i,nq + j) = d 2 Xk{npd,nqd)W ij > where Xfc is a 
continuous function and B' fc ) an. n X n matrix. We 
find by expansion about a c that there are n linearly 
independent matrices B^ with degenerate energies, of 
the form Bp 9 = e 2lxlvkl l n v q -\, where v is a fixed n- 
component vector for each I, n; e.g. for a c — 1/2 we have 



V = (^J I -l/y/2, y/l + 1/V2) 



for the minimum energy 

eigenvector. The subscript q — k wraps around mod n, 
and the Xk obey 
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Xk + V(x,y)xk = Exk, 

(4) 

where (l = (a — a c )w/md 2 and C depends only on l,n. 
This formula reduces to Eq. (0 for a c = 0/1. The com- 
parison in Fig. ^ shows excellent agreement of the wave 
functions for the cases a c — 1/2 and a c = 1/3. 

For this procedure to be consistent, the magnetic 
length £ on which Xk varies must be large compared to 
the small scale periodicity rid, but small enough that 
mTlyd ^C 1 for y ~ £. For a ID harmonic trap £1 = 
(m 2 fl 2 + rn 2 o; 2 /4C') _1 ' /4 and thus the analytic approxi- 
mation is valid for 



1 . 
— -r » [a - a c 



4C >(a 



a c ) 4 



(5) 



where (3 = md 2 Lu /ir is the dimensionless trap strength. In 
Fig. |2 we compare approximate analytical and exact nu- 
merical results for values of n < 8 as a function of a and 
[3. The plot shows large regions of validity of our calcu- 
lation around a c = (conventional FQH region) but also 
around a c = 1/2 and a c = 1/3. As expected for increas- 
ing n the regions of validity get narrower and are more 




1.5 



FIG. 2: Areas of validity of different fractions a c , in a 
one-dimensional harmonic trap. The solid (dotted) curves 
show regions with overlap between exact numerical wavefunc- 
tions and their analytical approximations larger than 99% 
(99.9%) and resulting energy difference smaller than 0.01/md 2 
(0.001 /rod ). The dashed lines are areas of validity obtained 
from Eq. taking ^> to mean a ratio larger than 50. 



easily destroyed by an external potential. The regions 
also get narrower for decreasing values of (3 with at the 
same time increasing number of n's fulfilling Eq. (J3J) at 
a = a c . Our model is thus only valid in a finite neighbor- 
hood of an a c if a suitable trapping potential is present. 
We extend the Hamiltonian to the many particle case 
by including an onsite interaction term into Eq. Q giving 
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ki,k 2 ,k3,k4 

where Xk(x, y) is now a bosonic field operator and the 
constants Gk u k 2 ,k 3 ,k 4 = Y.j w i-k 1 v ] -k 2 Vj-k 3 v ] - kA /n if 



ki+k2 = fc3+fc4 mod n and otherwise; e.g. for a c = 1/2, 
Gini ~ G2222 = 3/2 and G1212 = G2121 = G1122 = 
G2211 = 1/2. Particles with the same k interact more 
strongly since their v( fc ) are peaked on the same sites. 

For a c = 1/2 this effective Hamiltonian is equivalent 
to a bilayer FQH system with x± = Xi ± *X2 being 
the two "layers", and hence for weak potentials V the 
highest density zero interaction energy state is the 221 
state. This bosonic analogue of the fermion 331 state 
|l ll| with magnetic length £q = (mf2) -1 / 2 has filling fac- 
tor v = ptt£q = 2/3; it can be extended to larger n with 
v = n/(n + 1) by multiplication with (zi — Zj) 2 if par- 
ticles i,j are in the same "layer" and with (zi — Zj) if 
they are not. We expect a density step profile eventually 
breaking down into a vortex lattice, qualitatively simi- 
lar to the small a regime, with the 221 state being the 
lowest step. As long as Eq. (0 is a valid approximation, 
i.e. the particle spacing is much larger than nd, every 
step would occur at fixed v and not fixed v. From Fig. El 
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FIG. 3: Two-point functions: v = 1/2 Laughlin state (solid); 
v = 2/3 221 state gu = (722 (dashed) and g±2 = <?2i (dotted). 

the range of validity for a c = 1/2 is a — a c < 0.04 and 
hence for the v = 2/3 state, lattice filling pd 2 < 0.03. For 
U < J the energy gap is E g ~ Upd 2 which is about 10 
- 100Hz for typical optical lattice parameters consistent 
with [3 and will thus be measurable with near future 
technology. Measurement of the stepped density profile 
would confirm the existence of a series of incompressible 
states near a — a c of densities proportional to Q but 
would not confirm that they are FQH states. We will 
thus next investigate two methods feasible in optical lat- 
tices to identify and characterize these states further. 

We first study how shot noise measurements of the 
density correlations [lfl |ljj can reveal the two particle 
correlation functions g(|u>| 2 ) with w — x + iy and thus 
distinguish the FQH states from Mott insulating (MI) 
states for which such measurements were recently carried 
out [lg . As shown in [lfj results of these measurements 
are directly related to the expectation values 

n 

j2- R (fc;)- R *(fcj)- R *(fc;)- R (fcj) 



P2{i,j;i'j')k t k j -k i k j 



(7) 



where P2(hj',j' ' ,i')kik;kik- is the continuum 2-particlc 
density matrix. For LLL constant density states 
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and is thus determined by analytic continuation of the 



two-point correlation functions gk t k i\wi 



-Wj 



'). We have 



calculated g ki k ■ for the Laughlin and 221 states by Monte 
Carlo methods |lj|. The results in Fig. OUshow clear spa- 
tial antibunching up to distances of order £q which is very 
distinct from behavior of superfluid and MI states. 

Finally we calculate the Hall currents which can be 
used to characterize FQH states near a c . Unlike in a 
Galilean-invariant continuum FQH system the applica- 
tion of a linear potential V(x,y) = —may does not in- 
duce a particle current with velocity of exactly v x = 
a/20, since the lattice defines a rest frame. Instead, 




FIG. 4: Dimensionless Hall current I d — 1/(2 ' maVuQ/uj) 
against dimensionless number of atoms per unit length Nd — 
N/(2 3/2 m 3/2 m^/u~/uj) for Vi{y) = mui 2 y 2 /2, V 2 (y) = Vi(y)- 
may at small a and a, for no disorder (dotted straight line), 
maximal disorder (constant density of states) (black curve), 
and Lorentzian disorder of width umfl/10n (grey curve). 
Numbers under curve are filling factors at trap centre. The 
inset shows the v — 1/2 to v — 1 corner. 



when Eq. © is valid, the particle velocity is given by 
v x = 2Q,a/(AVt 2 + U! 2 /C) if the linear potential is super- 
imposed to a ID harmonic trap. The Hall current will 
thus cease for a = a c and for a < a c a counterintuitive 
negative Hall current flows. In a trapped system at equi- 
librium these currents will not be visible since they flow 
along equipotentials and steps in the density profile also 
lie along equipotentials. They can be made visible by 
putting the system out of equilibrium, e.g. by suddenly 
changing the trap from V\ to V2 as shown in Fig. 0] 

When sufficiently mild disorder is added to the optical 
lattice |2j, some of the particles become localized and 
cannot carry current, but the average velocity is still v x 
|llj . A simple model of this, valid for weak potentials, 
is to describe the disorder by a density of states Pd(E) 
giving the number of LLL states at energy E per unit 
energy interval per unit area and hence J pd{E)dE = 
ratt/ir. For a linear geometry we have 

M-Vi(y) 

J y J j> 

i"-Vi(y) 
N = L J &y J dEj2& ~ *3-i)Pd(E - Ej), (8) 



j>0 

with L the system's length, I the net current (both in 
the x-direction) and p, the chemical potential. The FQH 
states are at filling factors 9j and energies per particle 
Ej [for a w 0, Vj = (j + l)/2 and Ej w (umn/2n)j]. 

In a harmonic trap the disorder will not lead to FQH 
plateaus, only corners each time a new extended level be- 
gins to fill (c.f. Fig. [2J, but unlike the square-well case, 



it is possible to obtain the complete distribution Pd(E) 
by measuring 7" against N (or fi) and using Eqs. JSJ- 
The quasiparticles of our system carry fractional particle 
number J2tJ and fractional statistics [g] as in the conven- 
tional FQH system, but detecting this directly is likely 
to be more experimentally challenging. 

In summary we have investigated bosonic optical lat- 
tice FQH systems at large magnetic fields. We have 
worked out their ground states, magnetic lengths and 
the resulting fractional fillings, and found the Hall cur- 
rent to change its sign at a c . Furthermore, we have stud- 
ied ways for detecting these properties using recently de- 
veloped experimental methods. This work is supported 
by EPSRC through project EP/C51933/1 and QIP IRC 
(GR/S82176/01), and by the EU through OLAQUI. 
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